arXiv:1508.05786vl [cs.DS] 24 Aug 2015 


On the Displacement for Covering a (i— dimensional Cube 
with Randomly Placed Sensors^ 

Rafal Kapelko^* *’^’*, Evangelos Kranakis'’’^ 

^Department of Computer Science, Faculty of Fundamental Problems of Technology, Wroclaw University of 

Technology, Poland 

^School of Computer Science, Carleton University, Ottawa, ON, Canada 


Abstract 

Consider n sensors placed randomly and independently with the uniform distribution in 
a d—dimensional unit cube {d > 2). The sensors have identical sensing range equal to 
r, for some r > 0. We are interested in moving the sensors from their initial positions to 
new positions so as to ensure that the d— dimensional unit cube is completely covered, 
i.e., every point in the d— dimensional cube is within the range of a sensor. If the i-th 
sensor is displaced a distance di, what is a displacement of minimum cost? As cost 
measure for the displacement of the team of sensors we consider the a-total movement 
dehned as the sum Ma := some constant a > 0. We assume that r and n 

are chosen so as to allow full coverage of the d— dimensional unit cube and a > 0. 

The main contribution of the paper is to show the existence of a tradeoff between 
the d— dimensional cube, sensing radius and a-total movement. The main results can 
be summarized as follows for the case of the d— dimensional cube. 

1. If the d— dimensional cube sensing radius is 2 n^/d n = m‘^, for some m G 

N, then we present an algorithm that uses O total expected movement 

(see Algorithm]^ and Theorem]^. 

2. If the the a—dimensional cube sensing radius is greater than 2 n^ 

and n is a natural number then the total expected movement is O ^ ^ 

(see Algorithm]^ and Theorem]^. 

This sharp decline from ^) to 0(n^“ ^ (In n) ^) in the a-total movement of the 

sensors to attain complete coverage of the d— dimensional cube indicates the presence 
of an interesting threshold on the sensing radius in a d— dimensional cube as it increases 
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from to ( 3 i/d_i)( 3 i/d_i) ^;^- In addition, we simulate Algorithm|^and discuss 
the results of our simulations. 

Keywords: Displacement, Random, Sensors, d—dimensional Cube 


1. Introduction 

A key challenge in utilizing effectively a group of sensors is to make them form an 
interconnected structure with good communication characteristics. For example, one 
may want to establish a sensing and communication infrastructure for robust connec¬ 
tivity, surveillance, security, or even reconnaissance of an urban environment using a 
limited number of sensors. For a team of sensors initially placed in a geometric domain 
such a robust connectivity cannot be assured a priori e.g., due to geographic obstacles 
(inhibiting transmissions), harsh environmental conditions (affecting signals), sensor 
faults (due to misplacement), etc. In those cases it may be required that a group of sen¬ 
sors originally placed in a domain be displaced to new positions either by a centralized 
or distributed controller. The main question arising is what is the cost of displacement 
so as to move the sensors from their original positions to new positions so as to attain 
the desired communication characteristics? 

A typical sensor is able to sense a limited region usually defined by its sensing ra¬ 
dius, say r, and considered to be a circular domain (disc of radius r). To protect a larger 
region against intruders every point of the region must be within the sensing range of at 
least one of the sensors in the group. Moreover, by forming a communication network 
with these sensors one is able to transmit to the entire region any disturbance that may 
have occurred in any part of the region. This approach has been previously studied in 
several papers. It includes research on 1) area coverage in which one ensures moni¬ 
toring of an entire region nioiiia, and 2) on perimeter or barrier coverage whereby a 
region is protected by monitoring its perimeter thus sensing intrusions or withdrawals 
to/from the interior Cl E a i6i m. Note that barrier coverage is less expensive (in 
terms of number of sensors) than area coverage. Nevertheless, barrier coverage can be 
only used to monitor intruders to the area, as opposed to area coverage that can also 
protect the interior. 

1.1. Related work 

Assume that n sensors of identical range are all initially placed on a line. It was 
shown in 0 that there is an 0{n?) algorithm for minimizing the max displacement of a 
sensor while the optimization problem becomes NP-complete if there are two separate 
(non-overlapping) barriers on the line (cf. also ID for arbitrary sensor ranges). If the 
optimization cost is the sum of displacements then m shows that the problem is NP- 
complete when arbitrary sensor ranges are allowed, while an 0{'n?) algorithm is given 
when all sensing ranges are the same. Similarly, if one is interested in the number of 
sensors moved then the coverage problem is NP-complete when arbitrary sensor ranges 
are allowed, and an O(n^) algorithm is given when all sensing ranges are the same ifThl . 
Further, Q considers the algorithmic complexity of several natural generalizations 
of the barrier coverage problem with sensors of arbitrary ranges, including when the 
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initial positions of sensors are arbitrary points in the two-dimensional plane, as well as 
multiple barriers that are parallel or perpendicular to each other. 

An important setting in considerations for barrier coverage is when the sensors are 
placed at random on the barrier according to the uniform distribution. Clearly, when the 
sensor dispersal on the barrier is random then coverage depends on the sensor density 
and some authors have proposed using several rounds of random dispersal for complete 
barrier coverage ElIISl. Another approach is to have the sensors relocate from their 
initial position to a new position on the barrier so as to achieve complete coverage 
EJElKISl. Further, this relocation may be done in a centralized (cf. Em) or 
distributed manner (cf. il). 

Closely related to our work is 03, where algorithm MVi{n,y) was analysed. 
In this paper, n sensors were placed in the unit interval uniformly and independently 
at random and the cost of displacement was measured by the sum of the respective 
displacements of the individual sensors in the unit line segment [0,1]. Lets call the 
positions ^ for i = 1,2,... ,n, anchor positions. The sensors have the sensing 

radius r = ^ each. Notice that the only way to attain complete coverage is for the 
sensors to occupy the anchor positions. The following result was proved in El. 

Theorem 1 (cf. El l. Assume that, n mobile sensors are thrown uniformly and inde¬ 
pendently at random in the unit interval. The expected sum of displacements of all n 
sensors to move from their current location to the equidistant anchor locations ^ 
for i = 1,2,... ,n, respectively, is in Q(y/n). 


Algorithm 1 MVi {n, y) (Sensor displacement on a interval). 

Require: n mobile sensors with identical sensing radius r = ^ placed uniformly and 
independently at random on the interval [0, y]. 

Ensure: The hnal positions of the sensors are at the locations 
(so as to attain coverage of the interval [0, y].) 

1: sort the initial locations of sensors; the locations after sorting xi,X 2 , ■. ■ Xn, xi < 

X2< ■ ■ ■ < Xn- 

2: for * = 1 to n do 

3: move the sensor Si at position (^ — ^) 

4: end for 


In EIj Theorem[2was extended to when the cost of displacement is measured by 
the sum of the respective displacements raised to the power a > 0 of the respective 
sensors in the unit line segment [0,1]. The following result was proved. 

Theorem 2 (cf. El)- Assume that n mobile sensors are thrown uniformly and inde¬ 
pendently at random in the unit interval. The expected sum of displacements to a given 
power a of algorithm MVi{n, 1) is in 0 , when a is natural number, and in 

O when a > 0. 

An analysis similar to the one for the line segment was provided for the unit square 
in El- Our present paper focuses on the analysis of sensor displacement for a group 
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of sensors placed uniformly at random on the d—dimensional unit cube, thus also gen¬ 
eralizing the results of na from d = 2 to arbitrary dimension d > 2. In particular, our 
approach is the hrst to generalize the results of 1(131 to the d—dimensional unit cube 
using as cost metric the a-total movement, and also obtain sharper bounds for the case 
of the unit square. 

1.2. Preliminaries and notation 

Let d be a natural number. We dehne below the concept d-Dimensional Cube Sens¬ 
ing Radius which refers to a coverage area having the shape of a d-dimensional cubej^ 

Definition 3 (d-Dimensional Cube Sensing Radius). Consider a sensor 
Z{xi,x 2 ,...,xj,) located in position (xi, X 2 , ■ ■ ■, Xd), where 0 < xi,X 2 , ■. ■ ,Xd < 1- 
We define the range of the sensor area delimited by the d- 

dimensional cube with the 2^* vertices {xi ± r, X 2 ± r,... ± r), and call r the 

d-dimensional cube sensing radius of the sensor. 

We also dehne the cost measure a-total movement as follows. 

Definition 4 (a-total movement). Let a > 0 be a constant. Suppose the displacement 
of the i-th sensor is a distance di. The a-total movement is defined as the sum '.= 
Sr=i ^i- assume that, r and n are chosen so as to allow full coverage of the 
d-dimensional cube and a > O.J 

Motivation for using this cost metric arises from the fact that there might be a 
terrain with obstacles that obstruct the sensor movement from their initial to their hnal 
destinations. Therefore the a-total movement is a more realistic metric than the one 
previously considered for a = 1. 

In the analysis below we consider the Beta distribution. We say that a random 
variable concentrated on the interval [0,1] has the B(a, b) distribution with parameters 
a, b, if it has the probability density function 

B(a, o) 

where the Euler Beta function (see El) 

B(a, 6)= f x'^~^{l — x)^~^dx (2) 

do 

is dehned for all complex numbers a, b such as 3?(a) > 0 and 3?(6) > 0. Let us notice 
that for any integer numbers a,b> 0, we have 

B(a,6)-i = (3) 


^Recall that the generally accepted coverage area of a sensor is a d-dimensional ball. Our results can be 
easily converted to this model by describing a minimum d-dimensional ball outside this d-dimensional cube. 
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1.3. Results of the paper 

We consider n mobile sensors with identical d—dimensional cube sensing radius 
r placed independently at random with the uniform distribution in the d—dimensional 
unit cube (d > 2). We want to have the sensors move from their current location 
to positions that cover the d—dimensional cube in the sense that every point in the 
d—dimensional cube is within the range of at least one sensor. When a sensor is dis¬ 
placed on the d—dimensional cube a distance equal to d the cost of the displacement 
is d“ for some (fixed) power a > 0 of the distance d traveled. We assume that r and n 
are chosen so as to allow full coverage of the d—dimensional cube, i.e., every point of 
the region is within the range of at least one sensor. 

The main contribution of the paper in Section|^is to show the existence of a trade 
off between d—dimensional cube sensing radius and a-total movement that can be 
summarized as follows: 

1. For the case of the d—dimensional cube sensing radius ^n^/d n = m'^ for 
some m G N we present an algorithm that uses O (n^~ m ) total expected move¬ 
ment (see Algorithm]^ and Theorem]^. 

2. If the a—dimensional cube sensing radius is greater than 2 ni/^ 

and n is a natural number then the expected movement is O 
(see Algorithm]^ and Theorem]^. 

Notice that, for a = d Algorithm MVd{n, 1) uses 0{y/n) total expected movement 
while Algorithm LVd{n) uses 0(Vlnn) total expected movement. Therefore this 
sharp decrease from to 0{n^~^ (lnn)53) in the a-total movement of the 

sensors to attain complete coverage of the d—dimensional cube indicates the presence 
of an interesting threshold on the d—dimensional cube sensing radius when it increases 
from 1 to_ ^ _^ 

In Section plwe simulate Algorithmand provide the results of the simulations. 
Finally, Sectionfflis the conclusion. 

2. Displacement in d—dimensional cube 

Assume that n mobile sensors with the same d—dimensional cube sensing radius 
are thrown uniformly at random and independently in the d—dimensional unit cube 
[0, l]'^. Let a > 0 and d > 2. 

Our first result is an upper bound on the expected a—total movement for the case, 
where the d—dimensional cube sensing radius is 2 n^/d ■ Observe that in this case the 
only way for the sensors to attain complete coverage of the d—dimensional unit cube 
is to occupy the positions 

/ ^1 _ 1 ^2 _ 1 ^ \ 

^^1/d 2n^/d’ ^ nfld 2n^/dJ ’ 

for 1 < li,l 2 , ■ ■ ■ ,ld f: and li, I 2 ,... ,ld & Let us also notice that n = m‘^ 

for some m G N. We present a recursive algorithm MVd{n, 1) that uses O 
expected a—total movement. 
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Fig. 1: Nine mobile sensors located in the interrior of a unit square move to new positions according to steps 
(1 — 6) of Algorithm MV2(9,1). 


The algorithm is in two-phases. During the first phase (see steps (1 — 6)) we apply a 
greedy strategy and move all the sensors only according to the first coordinate. Figure 
[^illustrates the steps (1 — 6) of Algorithm MV 2 (n, 1) As a result of the first phase 
we have {d — 1)-dimensional cubes each with random sensors. Hence the 
first phase reduces the sensor movement on the unit d-dimensional cube to the sensor 
movement on the unit (d — 1)-dimensional cube. During the second phase (see steps 
(7 — 9)) we move sensors in the unit {d — 1)-dimensional cube. Notice that for the base 
case d = 1 we execute algorithm MVi{n, 1). 

We prove the following theorem. 

Theorem 5. Fix d G N. Let n = for some m G N and let a > 0. Assume that 
n sensors of d-dimensional cube sensing radius equal to o.re thrown randomly 

and uniformly and independently with the uniform distribution on a unit d-dimensional 
cube. The expected a—total movement of algorithm MVd{n, 1) is in O {v}~ ra ) . 

Proof We will prove the statement of the theorem by mathematical induction. Observe 
that the base case for d = 1 follows from Theorem|^[cf. ifT^ I. Let us assume the result 
holds for the number d— 1. Let a > 0. We will estimate the expected a—total movement 
at the steps (1 — 6). Let Xi be the fth order statistic, i.e., the position of the fth sensor 
in the interval [0,1] after sorting in step (1). It turns out (see |[T1) that Xi obeys the 
Beta distribution with parameters i^n — i -\- 1. We know that the density function for 
Xi (see Equation ([T])) is 


fxi {x) = i 





Therefore, the expected a—total movement in steps (1 — 6) of algorithm MVd{n, 1) is 
equal to 


E 


(a) 

(1-6) 


„i/d 

= E E 

j — l i—(j —+ l 




X 



1 

2n^/d 


x^-^{l-x)'^-^dx. 
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Algorithm 2 MVd{n, y) Sensor displacement on a d-dimensional cube when n is dth 
power of the natural number and d > 2. 

Require: n mobile sensors with identical d-dimensional cube sensing radius r = 
2 ^ 1 /d placed uniformly and independently at random on the d—dimensional cube 

[o,y]^ 


Ensure: The final positions of the 

(_yh _ y_ _yh _ yh yld _ 


sensors are at the locations 

V’-d ^ 1 < /i /n Id < 

2ndld ) I — ‘11 ''21 ■ ■ • I ‘d ^ 


and li,l2, ■ ■ ■ ,ld G N (so as to attain coverage of [0, y]‘^) 

1: Sort the initial locations of sensors according to the first coordinate; the locations 
after sorting = (a;i(l), 2:2(1),a;d(l)), 52 = (a:i(2), 2:2(2),2:^(2)),... 
Sn = (2:1 (n),2:2(71),.. .,Xd{n)), xi(l) < 2:1(2) < • • • < 2:i(n); 

2 : for j = 1 to do 
3: for 7 = 1 to do 

4: Move the sensor 5j_|_(j_i)„i/d at position 

( ^l/d 2 ndTd T ^i-\-(j—l)n^^d{ 2 ')^ . . . , 2:^_|_^j_i^^l/d (d)^ , 


5: end for 

6 : end for 

7: for j = 1 to do 

8: Execute MVd-i{n^‘^~^'>^^, y) for sensors 


9: end for 


Notice that, the expected a—total movement of algorithm MVi{n, 1) is equal to 

nl 




n' 

i / Jo 


i 1 




= E E 


2 = 1 


i 1 
'n“^ 


^(l-2:)"-*d2:. 


According to Theorem [^[cf. lfT2l l 


= O(ni-S), 


(4) 


when a > 0. 

Firstly, we estimate when a > 1. Notice that 

l 2 / + ^r<(|y| + kl)“<2“-i(|yr + |zr) for a>l,y,zGR. (5) 

This inequality is the consequence of the fact that f{x) = x°‘ is convex over i?+ for 
a > 1. Using Inequality § fory = 2^-(s-^)andz= (^- 2^)-(^-^) 
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we get 


X — 


< 2 ° 


J 


1 


2n^l'^ 


( 

X M 

a 

i j 1 1 

1 

1 

\n 2n J 


n nt/^ 2nt/‘^ 

2n 


( 6 ) 


We apply the definition of the Beta function (see Equation Q) with parameters i,n — 
i + 1, as well as Equation Q to deduce that 



x^-\i-xr-^dx 


i/d „(d-i)/d 


= E E 

j=l k=l 


1 


2n^/'^ 2n 


a n ' n 


lid ^(d-i)/d 


^E E 


/ 1 




v} d 


j=i fc=i 

Putting together Eormulas Q, (|^, and Q we obtain 

^( 1 - 6 ) ~ 0{n^~^) + 0(n^“3) = 0(n^“3), when a > 1, d > 2 
To estimate when 0 < a < 1 we define 


(7) 


( 8 ) 




j[ _ 

n^/d’ 2n^/'^ 


\l-x)^-^dx. 


Observe that. 


Fl 


^(i-e) “ E E " (*j)- 

j = l i=(j_l)„(rf-l)/£i + l 

Then, we use the discrete Holder inequality with parameters 7 and to derive 


,1/d 


E E 

7=1 d-1 

z=(j-l)n d +1 


(i.d) — 


d-1 
jn d 


E 

^ i=(d-l)r 





d-1 

i—[j — l)n d +1 


, V A \ 

E E (fS,)-| (")-“■ 


(9) 
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Next, we use Holder inequality for integrals with parameters ^ and and get 


X — 


X — 


J 


1 


j 1 


x''-\l-xY-\\ . da;< 


x^-\l-xr-H{ ]dx 


ni/'^ 2rdl^ ^ 

together Equation (|^ and Equation (|^ we obtain 


^{ 1 - 6 ) ~ ^ ‘^) > when 0 < a < 1. 


( 10 ) 


Observe that in step (8) of algorithm MVdin, 1) we have that mobile sensors 
are thrown uniformly and indendently at random in the unit (d — 1)-dimensional cube. 
According to inductive assumption the expected a—total movement at the step (8) is 
equal 0{{n~3 ~). Hence the expected a—total movement in steps (7 — 9) is 


in 




, 1 -* 


). Notice that the expected a—total movement 
in steps (1-6) is equal 0{'n}~^) (see Eormula ^ and Eormula (lOi). Therefore, the 
expected cost of displacement to power a of algorithm MVd{n, 1) is in 0{n^~^). 
This gives the claimed estimate for d and completes the proof of Theorem]^ □ 


Now we study a lower bound on the total displacement, when the d—dimensional 
cube sensing radius of the sensors is larger than ,^1/d ■ Eirst, we give a lemma which 
indicates how to scale the results of Theorem |^to d—dimensional cube of arbitrary 
length. The following lemma states that Algorithm MVd{n, y) uses O [y°‘ii}~ ra) ex¬ 
pected a—total movement. 

Lemma 6. Fix d G N. Let n = m'^ for some m G N and let a > 0. Assume that n 
sensors of d-dimensional cube sensing radius equal to 2 „^i/d thrown randomly and 
uniformly and independently with the uniform distribution on the [0, The expected 
a—total movement of algorithm MVd{n, y) is in O . 

Proof Assume that, n sensors are in the cube [0, yY- Then, multiply their coordinates 
by 1/y. From Theorem]^ the expected a—total movement in the unit cube [0,1]^^ is 
in O . Now by multiplying their coordinates by y we get the result in the 

statement of the lemma. □ 


A natural question to ask is: how to exploit the proposed Algorithm MVd{n, 1) 
when the number n of nodes is not a d—th power of natural number. Assume that n 
sensors have the d—dimensional cube sensing radius r = lyjrjd ■ To 

attain coverage of the cube [0,1]"^ choose sensors at random and use Algo¬ 

rithm MVd{n, 1) for the chosen sensors. Then similar arguments hold for Algorithm 
MVd{n,y). 

o3/d 

Notice that for / > ^- 3 i/d_x)( 3 i/d_i) we can do better. The following theorem 
states that Algorithm LVd{n) uses O I j expected a—total movement. 


9 


















Algorithm 3 LVd{n) Sensor displacement on a unit d-dimensional cube when d>2, 
p = j ( 2 + 2 ) , ^ = I (2 -I- 3 ) , xo is the real solution of the equation = 

3 such that xg > 3 

Require: n > [xqI mobile sensors with identical square sensing radius r = , 

o3/d 

f — (3i/d_i)(3i/d_ -jy placed uniformly and independently at random on the cube 
[ 0 , 1 ]'^. 

Ensure: The final positions of sensors to attain coverage of the cube [0, 

1 : Divide the d—dimensional unit cube into d—dimensional subcubes of side 


2 : if there is a d-dimensional subcube with fewer than 4 ^ ,^,d sensors then 
3 : choose sensors at random; 

4 : use Algorithm MVd{n, 1) that moves all n := sensors to equidistant 

points that are sufficient to cover the d-dimensional subcube; 

5 : else ^ 

6 : In each d-dimensional subcube choose (Alnn)^^‘^J sensors at random and 

use Algorithm MVd(n,t/) with n := (Alnn)^^‘^| ,y '■= 1 , p/d\ to move 

L -I J 

the chosen sensors to equidistant positions so as to cover the d-dimensional sub¬ 
cube; 

7 : end if 


sensors at random and 


10 



Theorem 7. Fix d G TV \ {1} and a > 0. Let f > and n > [xol, 

where Xq is the real solution of the equation 9 ^ 2 +s) in a “ ^ such that xq > 3. Assume 

that n sensors of d-dimensional cube sensing radius r = 2 rS/d thrown randomly 
and uniformly and independently with the uniform distribution on the [0, l]"^. The ex¬ 
pected a—total movement of algorithm LVd{n) is in O ( 


Proof Assume that d G TV \ {1} and a > 0. Let p = | (2 + g) and A = | (2 + |) , 
Xq is the real solution of the equation 9 , . — = 3 such that xq > 3. First of all, 

observe that ^ j > 3 for n > [xq]. We will prove that Algorithm LVd(n) uses 

O ) expected a—total movement. There are two cases to consider. 

Case 1: There exists a d—dimensional subcube with fewer than 

1 n 

3 ' 


( " ^ 
\p\nn I 


1 /d 


sensors 


. In this case choose \ sensors uniformly and randomly from n sensors. 


Applying the inequalities [xj > x - 1 and / > ( 3 i/d_i)( 3 i/d_i) > J/dli we deduce 
that 




,1/d 


> 


f/d _ I 31/d 
n^/d 31/d _ I 


> 1 for n > 3. 


Therefore, the \ chosen sensors are enough to attain the coverage. The expected 

a—total movement is O = O (n^“^) by Theoremj^ 


Case 2: All d—dimensional subcubes contain at least 4 - 


the inequality [xj < x we deduce that, 

(Alnn)^^‘^J 




sensors. From 


d 1 
< - 
- 3 


( " ^ 
\ p In n J 


1 /d 


d ' 


Hence it is possible to choose 


(Alnn)^^‘^J 


sensors at random in each d-dimensional 


subcube with more than 4 ^- " ,,., d sensors. Let us consider the sequence 


33/d 


an = 




1 


,1/d 


( 3 l/d-l)( 3 l/d-l) L 
for n > [xq] . Applying inequality [xj > x — 1 we see that 

/... , 1 /d 




l/d 


dn ^ 


(3l/d-l)(3l/d_l) 


Alnn)"/ -1 


) , 2 , 1 /d ^ 


plnn^ 


p\nn 


l/d 


- 1 
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32 /d 


(3i/d-l)(3i/d_i) 


1 - 


(A Inn 


,1/d 


1 - 


plnn\ 


Observe that 


plnnl 1 __ 

- < - < o for n > [xol 

n 6 A in n o 


Putting together Equation (111 and Equation ([T^ we get 


( 11 ) 


( 12 ) 


(Alnn)^^^J 


i/di'^Z" 




1 /d 


p\nn 




> a“ > 1. 


Therefore, (Alnn)^^‘^J chosen sensors are enough to attain the coverage. By the 

independence of the sensors positions, the (A In n) chosen sensors in any given 
d-dimensional subcube are distributed randomly and independently with uniform dis¬ 
tribution over the d-dimensional subcube of side y = j -— ^ ^ ^ gy Lemma the 

expected a—total movement inside each d-dimensional subcube is 


O 


I 


\ 


( " ^ 
\ p In n J 


1 /d 


( 2 llnn)r/'^J' 


\ 


= O 




{In n ): 


nd 


-(Inn) ) . 


Since, there are 


( ri \ 
\ p In n J 


1 /d 


d-dimensional subcubes, the expected a—total move- 

It remains to 


ment over all d-dimensional subcubes must be in O {n^ 2 d (ly^) . 

consider the probability with which each of these cases occurs. The proof of the theo 
rem will be a consequence of the following Claim. 

Claim 8. Let p = | (2 -I- ()) . The probability that fewer than 






^ sensors 


fall in any d-dimensional subcube is < 

II 

Proof. (Claim|^ Eirst of all, from the inequality [a;J < a; we get 




( 2 + 5 ) Inn 

( " 

n 

\p\nn J 


2 

< -. 
- 3 


Hence, 


d< 


\plnny \p\nnj A \plnny 


(2-1- l) nlnn 


(13) 
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The number of sensors falling in a d-dimensional subcube is a Bernoulli process with 
probability of success ^ ^ . By Chernoff bounds, the probability that a given 

d-dimensional subcube has fewer than 


\plnnj \ \ plnn J 


(2 + 2 ) nlnn 


sensors is less than e 


-(i+m) In" < -J- 


-. Specifically we use the Chernoff bound 


Pr[X < (1 - 6 )m] < 


m = 


L(t 


/ ( 2 + 3 ) In" 



V p In n } 


d 


. As there are 


( " 1 

\plnn J 


1/d 


d-dimensional subcubes, the event that one has fewer than 


( " ^ 
\ p In n J 


1/d 


\ 


(2 + 2 ) nlnn 


( ^ 

\ p In n J 

K^n 


l/d 


d ' 


sensors occurs with probability less than —j_ Equation ( [T3] l com¬ 
pletes the proof of Claim|^ □ □ 

Using Claimj^we can upper bound the expected a—total movement as follows: 


/ 


1 - 


( n \ 
\ p In n J 


1/d 




n^+2d 


( 


O 2d 


Inn 

n 


( ^ 
\ p In n J 


1/d 




2 d 


O 2d ^ = 




O V} 2d 


Inn 

n 


which proves Theorem]^ 


□ 


□ 


3. Simulation Results 

In this section we use simulation results to analyze how random placement of sen¬ 
sors on the square impacts the expected a—total movement. 

We repeated 3 times the following experiments. Firstly, for each number of sensors 
n G {2^, 3^, 4^,..., 60^} we generated 32 random placements. Then we calculated 
the expected a—total movement according to Algorithm MV 2 {n, x). Let £^n ,32 be the 
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Fig. 2: The expected a—total movement of Algorithm MV 2 {‘n, 1). 


average of 32 measurements of the expected a—total movement. Then, we placed the 
points in the set {(n, En^ 32 ) : n = 2 ^, 3^, 4^,..., 60^} into the picture. 

Figure [^illustrates the described experiments for Algorithm MV 2 {ri, 1) when a = 
2 and a — 4. The additional line in the above pictures is the plot of the function 
which is the theoretical estimation. Black dots which represent numerical results are 
situated near the theoretical line. According to the proof of Theoremj^the steps (7-9) of 
Algorithm MV 2 (n, 1) conctribute the asymptotics. Notice that, the expected a—total 
movement in steps (7-9) of Algorithm MV 2 {n, 1) is equal to 


E 


(a) 

(7-9) 







Applying the Formulas for and E^p_g^ in any mathematical software that per¬ 

forms symbolic calculation we get 


E 


( 2 ) 

(7-9) 


1 

6 


y/n and E 


(4) 

(7-9) 


1 

To' 


Therefore, M 2 = ^y/n and M 4 = E. 

Figure ^illustrates the described experiments for Algorithm MV 2 (n, |) when a = 
2 and MV^n, 2) when a = 4. The additional line in the above pictures is the plot of 
the function which is the theoretical estimation. Black dots which represent numerical 
results are situated near the theoretical line. According to the proof of Lemma [^ we 
have M 2 = (lY \y/n = j^y/n and M 4 = 2 ‘^E = 4^ 


4. Conclusion 

In this paper we studied the movement of n sensors with identical square sensing 
radius in d dimensions when the cost of movement of sensor is proportional to some 
(fixed) power a > 0 of the distance traveled. We obtained bounds on the movement 
depending on the range of sensors. 
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Fig. 3: The expected a—total movement of Algorithm MV 2 {n, y). 
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